Aerospace  Environmental  Support  Center 
Technical  Memorandum  70-3 
IONOSPHERIC  ELECTRON  DENSITY  PROFILE  MODEL 

July  1970 

by 

Thomas  D.  Damon 
and 


Franklin  R.  Hartranft 


BLANK  PAGES 
IN  THIS 
DOCUMENT 
WERE  NOT 
FILMED 


TABLE  OF  CONTENTS 

Fa  Re 

I.  Preface .  1 

II .  Introduction . . . .  2 

III.  Development  of  the  Model  . .  3 

IV.  Description  of  the  Computer  Program . . .  ^ 

V.  Description  of  the  Computer-Produced  Profile  . .  7 

VI.  Evaluation  of  the  Model  . . 7 

VII.  Summary  and  Conclusions  .  9 

VIII.  References  . 10 

Figures  1  -  9  . . . . .  11 

Appendix  A.  Computer  Program  ''MODEL"  .  21 

Appendix  B.  Sample  Compute*  Output  .  31 


4 


i.  Preface. 


Recent  investigation  ol  the  effects  of  ionospheric  retardation  and  refraction  on 
satellite  tracking  radars  has  generated  a  need  for  a  means  to  predict  the  errors  and 
correct  for  them.  This  paoer  describes  e  project  undertaken  by  ilth  Weather  Wing  to 
produce  a  realistic  electron  density  proflic  based  upon  parameters  which  can  be 
forecast  reasonably  accurately.  Tne  authors  wish  to  acknowledge  the  help  provided 
them  ir  this  project.  Lt  Colonel  Hansrote  provided  tne  impetus  for  producing  such 
a  model.  Capt  Jack  Wro'oel  solved  our  initial  problems  of  scale  height  by  providing 
"Wrobel's  Equation."  MSgt  Birch  and  TSg*.  'uster  analyzed  and  evaluated  the  model 
against  actual  observations.  Mrs.  Green  accomplished  the  manuscript  typing.  Thanks, 
also  to  Lt  Bo  Eross  for  his  system  analysis  suggestions. 


Thomas  D.  Damon,  Major,  USAF 
Technical  Development  Branch 
Detachment  1,  4  WWg 


Franklin  R.  Hartranft,  Captain,  USAF 
Autonation  Branch 
Aerospace  Sciences  Division 
Hq  WWg 


II. 


Introduction. 


'.’he  development  of  a  computer  program  for  predicting  electron  density  profiles 
was  prompted  by  the  realization  that  ionospheric  retardation  and  refraction  pro¬ 
duced  errors  in  range  and  azimuth  of  satellite  tracking  radars.  These  errors  are 
of  the  same  order  of  magnitude  as  those  produced  by  tropospheric  effeccs  when  the 
UiiF  rauars  are  operating  above  a  few  degrees  elevation.  Since  the  effects  are 
rather  small,  it  was  assumed  that  a  simple  model  from  100  km  to  1000  km  would  be 
sufficient.  However,  as  development  work  began,  other  requirements  for  electron 
density  profiles  became  apparent.  A  three-dimensional  ionosphere  for  HF  ray  tracing 
wi ict.  requires  considerably  more  accuracy  in  the  lower  ionosphere,  was  requested. 

In  addition,  total  electron  content  for  correcting  for  Faraday  rotation  in  some 
navlgat ioral  satellites  requires  a  model  extending  higher  than  1000  km. 

The  program  described  in  this  report  has  been  used  routinely  for  about  eight 
months  to  predict  electron  density  profiles  for  the  FPS-85  radar  at  Eglin  AFB,  Fla. 
Results  are  encouraging  enough  to  warrant  publication.  It  should  be  considered  an 
interim  report,  however,  as  improvements  are  sure  to  be  required  as  its  accuracy 
is  evaluated  for  different  purposes. 


? 


III.  Development. 

The  ionospheric  electron  density  profile  model  presented  in  this  paper 
consists  of  the  sum  of  three  Chapman  layers  (E,  FI,  and  F2).  Each  layer  is  of 
the  form 


Nh  -  Nmax  exp  a[l-Z-exp(-Z)] 
where  Z  -  0"‘-hmax)/hg 

*  electron  density  at  height  h 

Nmax  “  electron  density  at  the  peak  of  the  Chapman  layer  hmaX 

ha  «  scale  height  at  the  peak  (except  for  the  topside  of  the  F2  region) 

The  value  of  the  constant,  a,  depends  upon  whether  electrons  are  lost  by 
attachment  or  by  recombination.  While  neither  process  is  unique  in  any  layer,  a 
is  assumed  to  be  0.5  for  the  E-layer  and  1.0  for  the  FI  and  F2  layers. 

Electron  densities  in  the  topside  ionosphere  are  controlled  by  complex 
motions  rather  than  a  production-loss  balance  and  cannot  be  successfully  described 
strictly  by  a  Chapman  layer.  An  effort  was  made  to  keep  from  over-complicating 
the  model  and  still  obtain  the  best  topside  profllq.  After  some  experimentation 
a  fit  was  obtained  by  simply  using  the  Chapman  equation  for  the  topside  ionosphere, 
but  computing  the  electron  densities  by  using  a  variable  scale  height  throughout 
the  region. 


The  scale  height  profile  is  calculated  from  the  equation 


log  h 


2.186  X  lO-2 


203.  ^7 


This  equation  describes  the  scale  height  of  a  simple  standard  atmosphere  and  was 
derived  by  Capt  J.  Wrobel  (private  communication). 

Critical  frequencies  for  the  E  and  FI  regions  are  determined  from  regression 
equations  [ 1 3 ,  [2]. 

f0E  ■  0.9CU80  +  1 .  MUR)  cos  xl0,25 

f 0F1  -  1.26f0E  +  0.5 

where  R  *  the  twelve-month  running  mean  sunspot  number 
x  ■  the  solar  zenith  angle 

When  x  exceeds  90°,  fQE  is  set  to  0.7  MHz.  When  x  exceeds  135° »  f0E  is 
set  to  0.3  MHz. 

The  F2  region  critical  frequency  may  be  predicted  from  the  ITS  (ESSA) 

coefficients  by  predicting  a  sunspot  number  (R)  [1].  It  may  also  be  predicted 

manually  on  a  short-term  basis  by  the  Air  Force  Aerospace  Environmental  Support 
Center.  For  post  analysis  purposes,  an  observed  value  may  be  used. 

The  height  of  the  peak  of  the  E  region  is  assumed  to  be  120  km.  After  some 

experimentation,  the  FI  peak  was  placed  halfway  between  the  E  and  F2  peaks. 

The  height  of  the  f2  peak  is  calculated  by  using  Shimazaki's  equation  [3]: 

11)90 


where  the  M(3000)  factor,  M,  may  be  predicted  in  a  manner  si'..ii]ar  tc  the  prediction 
of  foF2,  or  observed.  Computations  of  hmax  using  M(3000)  were  found  to  be  accurate 
within  20  km  at  mid  latitudes.  If  a  more  accurate  measure  of  hmax  is  available, 
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such  as  hp  ^2,  an  artificial  M(3C00)  may  be  calculated  from  the  Shimazaki  equation 
and  used  as  an  input  into  the  computer  program. 

IV.  Description  of  the  Computer  Program. 

A  copy  of  the  computer  program  used  to  compute  an  electron  density  profile  is 
listed  iii  Appendix  A.  The  program  is  written  in  IBM  7090  FORTRAN  IV.  There  are 
three  input  options  (all  of  which  are  concerned  with  the  method  of  obtaining  foF2 
and  M( 3000;).  Two  output  options  are  available,  depending  upon  the  representation 
of  the  profile  required. 

The  program  computes  electron  densities  independently  for  each  of  three  regions 
(E,  FI  and  F2).  The  base  of  the  profile  is  J 00  km  and  computations  are  made  at 
5  km  increments  to  1000  km.  The  three  regions  are  added  toge;  jr  to  give  the  total 
electron  density  at  each  increment  of  altitude.  Electron  density  is  not  permitted 
to  decrease  with  altitude,  but  is  held  constant  across  "valleys"  in  the  profile. 

Total  electron  content  in  a  one  square  meter  cross  section  up  to  a  given 
altitude  is  also  computed.  An  initial  electron  content  is  established  at  95  km,  to 
represent  the  total  content  below  100  km.  A  calculation  of  plasma  frequency  is  made 
from  the  electron  density  f^r  each  5  km  interval. 

Input  parameters  are  read  from  data  cards.  The  first  data  card  indicates  out¬ 
put  options  (Table  1).  The  second  data  card  contains  information  pertaining  to  the 
geographic  location  of  the  profile.  The  format  of  the  card  is  the  same  for  all 
input  options  (Table  2). 

As  previously  mentioned,  there  are  three  input  options  which  determine  the 
method  by  which  foF2  and  M(3000)  are  introduced  into  the  program.  Data  card  number 
3  contains  information  pertinent  to  the  profile,  including  the  input  option  vari¬ 
able  IOFT.  Table  3  lists  the  input  parameters  on  the  third  data  card  and  indicates 
which  of  them  are  used  by  the  program  under  each  of  the  three  ■'nput  options.  If 
IOPT  *=  1,  foF2  and  M(3000)  are  computed  from  a  card  deck  of  ITS  Prediction  Coef¬ 
ficients.  (Subroutines  used  to  compute  foP2  and  M(3000)  from  ITS  coefficients  were 
extracted  from  a  program  published  in  [4].)  If  IOPT  ■  2,  a  long-term  data  tape 
containing  sunspot  dependent  coefficients  of  foF2  and  M(3000)  is  read  to  determine 
foF2  and  M(3000).  Finally,  if  IOPT  ■  3,  foF2  and  M(300Q)  are  read  explicitly  from 
the  data  card. 
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TABLE  1 


Card  Column 

Variable 

Explanation 

1 

I  PLOT 

If  I PLOT  -  1,  a  profile  of 
plasma  frequency  vs  height  Is 
plotted.  If  IPLOT  *  0,  plot 
is  suppressed. 

2 

IPNCH 

If  IPNCH  -  1,  the  17  most 
significant  points  depicting 
the  profile  are  punched  onto 
data  cards.  If  IPNCH  >*  0,  the 
punch  routine  is  suppressed. 

TABLE  2 

Card  Column 

Variable 

Explanation 

1-6 

CLAT 

Latitude 

7 

NORS 

Hemisphere  (N,  S) 

8-13 

CLONG 

Longitude 

lit 

IHEM 

Hemisphere  (E,  W) 

15-38 

NAME 

Name  of  Station 

5 


» 


TABLE  3 


Card  Column 

Variable 

Required  Under  Option 

Explanation 

1-2 

IYR 

1.  2,  3 

Year 

3-4 

MNTH1 

1.  2,  3 

Month 

5-6 

MNTH2 

1,  2 

Used  if  mean  of 
two  months  coef¬ 
ficients  are 
required. 

OO 

IDA 

3 

Day 

9-12 

IBHR 

1,  2,  3 

Beginning  time 
of  set  of  con¬ 
secutive  pro¬ 
files  (GMT). 

13-16 

IEHR 

1,  2,  3 

Ending  time  of 
set  of  consecu¬ 
tive  profiles 
(GMT). 

17-18 

INC 

1,  2,  3 

Increment  of  time 
step  (hours). 

19-21 

JDAY 

1,  2,  3 

Julian  Day 

22 

IOPT 

1,  2,  3 

Option 

26-30 

SSN 

1,  2,  3 

Sunspot  Number 

31-40 

F0F2 

3 

foF2  (explicit ) 
MHz  and  tenths 

41-50 

EM300Q 

3 

M( 3000)  Explicit ) 
hundredths 

51-56 

IVB 

1,  2 

Beginning  of 
valid  time  (i.e.  , 
10  May). 

57-62 

IVE 

1,  2 

Ending  of  valid 

time  (i.e.,  20 
May ) . 


NOTE:  All  numbers  are  integers  except  SSN,  FOF2,  and  EMjOOO.  These  three  are 
floating  point,  punched  with  a  decimal  point,  anywhere  in  the  field. 


Y 
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V.  Description  of  the  Computer  Produced  Profile. 

Appendix  2  Is  a  sample  profile  produced  by  computer.  The  profile  Is  In  four 
pectlons.  The  first  section  provides  a  summary  of  Input  da. a  and  pertinent  infor¬ 
mation  for  each  of  the  three  regions.  The  second  section  Is  the  profile  Itself, 
listing  values  of  height,  E-reglon  density,  Fl-reglon  density,  F2-region  density, 
total  density,  cumulative  electron  content,  plasma  frequency  and  scale  height  for 
each  5  km  Increment  of  the  model.  Output  of  the  third  section  depends  upon  the 
value  of  the  output  option  IPLOT  (see  Table  1),  and  plots  a  graph  of  plasma  fre¬ 
quency  vs  height  f>r  the  model.  The  fourth  and  final  section  depends  upon  the  value 
of  the  output  option  IPNCH.  If  selected,  the  17  most  significant  values  of  plasma 
frequency  describing  the  profile  are  chosen  objectively  and  written  onto  magnetic 
tape  for  punching  onto  data  cards.  In  addition,  a  checklist  of  the  points  selected 
Is  printed. 

VI.  Evaluation  of  the  Model. 

An  evaluation  of  this  model  was  made  by  comparing  with  observed  electron 
density  profiles  and  with  total  electron  content  measurements. 

Figures  1  to  8  show  model  monthly  median  profiles  for  Wallops  Island,  Va. , 
during  1968,  compared  with  the  observed  profiles  available  from  World  Data  Center 
A,  Boulder,  Colorado.  Excellent  results  are  obtained  during  winter  and  at  night. 

The  July  loOOZ  (mid-day)  Is  the  worst  case  among  severa^  dozen  such  comparisons  at 
various  locations  and  times. 

In  Table  4,  the  total  electron  content  calculated  to  1000  km  compared  with 
observations  of  total  content  from  Bedford,  Mass.,  to  geostationary  atellltes,  a 
path  which  passes  through  the  F2  peak  near  Wallops  Island.  These  observations, 
courtesy  of  Jack  Klobuchar,  Air  Force  Cambridge  Research  Laboratories,  are  converted 
to  vertical  Incidence  by  assuming  a  cosine  correction  factor.  As  expected,  the 
model  Is  generally  lower  than  the  observations  since  It  cuts  off  at  1000  km.  It  Is 
Interesting  to  note  that  in  the  summer  daytime,  when  the  model  overestimates  the 
bottomslde  content  (Figure  6),  it  underestimates  the  total  content.  This  Implies 
that  more  electrons  are  present  in  the  topside  than  the  model  predicts. 

A  third  comparison  Is  shown  In  Figure  9.  Here,  total  content  to  1000  km  from 
the  model  Is  compared  with  the  total  content  on  near  vertical  Incidence  paths  to 
synchronous  satellites  In  the  vicinity  of  Hawaii  [5].  The  shape  of  the  diurnal 
curve  Is  good  and  the  results  are  again  excellent  at  night  but  are  underestimated 
at  midday. 
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TABLE  4 


TOTAL  ELECTRON  CONTENT 
<10-  M' c ) 
WALLOPS  ISLAND  1968 


GMT 

OBSERVED  TEC 

MODEL  TO  1000  KM 

PERCENT  DIFFERENCE 

January 

0100 

1.2 

1.0 

-10 

0500 

.66 

.56 

-15 

2000 

4.1 

3.5 

-15 

March 

0200 

1.6 

1.4 

-13 

0500 

.98 

.95 

-3 

1100 

.65 

.61 

-6 

1600 

3.3 

3.2 

-3 

2100 

3.7 

3.^ 

-8 

May 

0300 

1.5 

1.4 

-7 

0700 

.85 

.93 

+9 

1200 

1.4 

1.2 

-14 

1600 

2.3 

2.0 

-13 

2200 

2.7 

2.1 

-22 

July 

0400 

1.2 

1.1 

-8 

0800 

.62 

.60 

-3 

1300 

1.4 

1.3 

-7 

1800 

1.8 

1.6 

-11 

2300 

2.1 

1.6 

-24 

September 

0100 

1.4 

1.4 

0 

0500 

.91 

.89 

-2 

1000 

.42 

.46 

+  10 

1500 

2.4 

2.4 

0 

2000 

2.9 

2.7 

-7 

November 

0200 

.75 

.65 

-13 

0600 

.55 

.52 

1100 

.39 

.36 

-8 

1600 

3.2 

3.3 

+3 

2100 

3.2 

2.8 

-13 
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VII,  Summary  and  Conclusion. 

The  electron  density  profile  model  came  about  as  an  attempt  to  produce  a 
reasonably  simple  method  of  predicting  electron  densities  In  the  100-1000  km  range. 
This  model  should  not  be  considered  finul,  by  any  means,  even  for  the  current 
applications.  The  scale  height  profile  should  be  Improved  to  Include  diurnal  and 
seasonal  variations.  The  model  should  be  extended  from  Its  upper  limit  of  1000  km 
to  the  plasmapause.  Variations  ol  electron  density  with  geomagnetic  activity 
should  be  included.  Improved  prediction  or  specification  of  any  of  the  Input 
parameters  will,  of  course,  Improve  the  accuracy  of  the  profile  model.  The  height 
of  the  F2  maximum  Is  probably  the  most  Important  of  these. 
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Computer  Program 
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21 


UBFTC  MQbfcL 

C IMENS ION  SCALE! 182). IHGT! 182) 

_  DIMENSION  VAH91) .NAPE (4) ,STOK( 182) 

"COMMON  78/  K  ( 1 A »  ,U  (17,76 ) ,  KX  (  M I ,  UX  ( 1 7 . 76 )  fKlUA)  »KX1  ( 14)  , 
IU1(17.76).UX1(17.76) 

UATA  "IE,  IS/lHfc,  1HS/ 

C 

SET  UP  CONSTANTS  FOR  KAO  IAN  CONVERSION  #«•••*• 

_C 

"c 180*3. 14 i 592 f 
P I  2*0180/2.0 
AK*C 18C/ 180 .0 
UK*lto0.0/Cl80 
GLT* 1.36135662 

__  GL_G=  1.22 1.730  30  _  _ _ 

KARL=G 

C 

SCALE  HEIGHT  COMPUTATION  (WROBEL-S  EQUATION)  '  ***** 

C 

H*V5T  ~ 

CO  10  I_f  1.182. _ _  _ 

SCALEU  »=TAL0G(~H)/2.186E-b2f-2G3.AA7~ 

IHGT1I )*H 

1 F ( I . EG .1 )  SCALE!  I  >  =  6.6 


10  H=H+5.0 
C 

C***»**P EAO  OUTPUT  OPTIONS  ****** 

c  .  . 

REAUJJ>,  m  I  PLOT ,  I PNCH 

11  . FORMAT  ("2 ID”  '  . 

C 

C*i*i**REAO  LOCATION  OF"  STATION  ****** 

C 

99V  RLAU  (  5 . aT  CL AT  ,  NOR  S  ,"C LONG  . THEM , !  NAME  fl ) »  1*1  ,A) 
A  FORMAT  (2(F6.1«A1),AA6) 

c  “ . .  .  ~  '  . . 

C******  ChECK  FOR  BLANK  CARU  (END  OF  ALL  DATA)  *••• 
C  ‘ 

I F ( CLAT • EQ.0.0 )  GO  TO  777 

c . . . .  . . . .  "  “  . 

C******  CHtCK  FOR  EAST  CRN  OR  SOUTHERN  HEMISPHERE  *« 
C 

IF(NORS.EQ.IS)  CL AT*- 1. *CLAT 

IF! IhE¥VEQ.IE)  CLONG*- 1.C«CL0NG  " 


8000  REAU  (5.AA)  I YR , MNTH1 , MNTH2, ILA , I BHR, IEHR, INC . JDAY . I  OPT , SSN. FCF 
ITfOGOG.TvBTrvE  "  - - - - ~ - 


AA  FORMAT  ( A 1 2 , 2  I A, 1 2, 1  3 , 1 1 , 3X, F 5 . 1 , 2F 10. 5 , 2A6 ) 
iHR* IBHR  " 

INC* INC* 100 

c . 

C******  CHECK  OPTION  1  *  ITS  COEFFICIENTS  READ  FROM  CARUS 

IT"  . .  . . 2"”*  long- Term  "o’4Ta"Tapt  . . 

C  3  =  FCF2  AND  MI30C0)  EXPLICIT  *••** 

C 


IF( IOPT.EQ.l. AND. KARO. fc'Q.C)  GO  TO  1 
I F ( I  OPT •  EQ • 1 •  ANU.KARD.fcO. 1 )  GO  TO  34 
JFJ  ICPJj.EQ.2i  GO  TO  2 

lpi  lOPT.Ec.3J  'GO  'TO  ~3  ""  '  "  " . 

GO  TO  33 
C 

. . *  I  OPT  *  i,  READ  CAROS  . 

C 

1  CALL  KEACUU'U) 

CALL  READUUX.UXV 
I F ( MNTH2.GT .0 )  CALL  READU(K1,U1) 
IF(MNTH2.GT.0)  CALL  REAOU(KXliUXl)  ' 
KARu* 1 
GO  TO  34 
C 


£•*****”  10PT  =  2»  READ  LONG-TERM  DATA  TAPS'  ’ •••••• 

C 

2  CALL  L T APE  ( MNT H l » S  SN ,  K ,  U ,  KX , UX ) 

1 F  (  MNTH2 • GT.O )  CALL  LT  APE  (  MNTH2 , S SP,K1  ,«U1 , KX1  ,UX1 ) 
34  CONTINUE' 

800  CALL  001T(CLAT,CL0NG, IHR , F0F2 , EM3G00 , MNTHl ,MNTH2 ) 

GOTO  f  '  '  . . . 

33  WRITE  (6,333) 

333  Format  ! Thi, I3HERR0K  ;n  Topri 

GO  TO  777 
J'  CONT I NUt 

C _  _  _  _  _ 

'~CY. ”iOPf  =  3,  BEGIN  COMPUTATIONS  ' 

C 

C 

C*«*«*«  NOW  HAVE  F0F2  AND  M30C0  BY  ONE  OF  THREE  METHODS 

c  . 

C******  CALCULATE  SOLAR  ZENITH  ANGLE  «•«•«« 


IFRST=IHK/10C 
SECND* IHR- I FRST* ICC 
SECNU=SECND/60. 

GMT* FLOAT ( I FRST ) ♦SECNC 
IF( GMT. EC-0.0)  GMl *24 • 

0= JDAY 

SSP*-23.45*C0$( (U+ 10. )/365.»C 180*2  .J 
SSP=SSP*AK 
SSL=I5.0«GMT-180.0 
Z  = ( SSL-CLONG  j  *AK 

COMP*  S I  N(CLAT*AK)«SlN(SSP)+COS(CLAT*AK)*CQSCSSP)*COS(Z) 

- caMWR'co'siccw)  ■  - - - 

C0MP*A8S(C0MP) 

RANG=COMP 

ZANG=C0MP*BK+0.5 

C 

C  CALCULATE  GEOMAGNETIC  LATITUDE 


24 


C 

C« 

C 

c 

C_. 

c 

c 

c» 

c" 

c" 

c» 

c 


c« 

c 

c 

C. 

c" 

c 

c* 

c 


c 

c« 

c' 

c 

c* 

c 

c 

c« 

c 

c 

c« 

c 

c 

t« 

c 


GAT=ARCO$ ( S1N(GLT)»SIN(CLAT*AK)+C0SIGLT)*CUS(CLAT*AK) 
1*CCSICL0NG*AK-GLG) ) 

GL  ATMPI2-GAT  )*WK 
R*  S5N  . 

•  '  NC'h  have  needed  parameters  for  equations  **•*•* 

_ HEIGH r  OF_E-REGl ON  SET  TO  120  KM  * * 

HE* 120.0 

••  COMPUTE  SCALE  HEIGHT  OF  E-REGICN  •••••* 

TE  =  (  AL0G(HE)/2.18St-02)-2n';.447  _ _ 

**  COMPUTE  FOE  *»•*** 

PAR  T  *0. ( 160.0*l.44*R)*C0SIRANG) ) 

I  F ( PAR  T  .  GE . 0 . )  FC  E *  P  AR  T • • 0 • 2  5 
I F (PART «L  T .0. )  FQE*0. 7 

*•  IF  SOLAR  ZENITH  ANGLE  GREATER  THAN  130  DEG »  FOE  SET  10  0 
IF ( Z ANG .GE.90.0 )  F0E*G.7 

IF  SCLAK  ZENITH  ANGLE  GREATER  THAN  90.  DEG  t  FOE  SET  TO  0. 
IF( ZANG.GC. 130. )  FOE  =  C. 3 
**  COMPUTE  F0F1 

F0Fl»1.26*F0E+0.5  _ 

**  COMPUTE  MAX  DENSITY  OF  E-REGION  •***«* 
ENE*1.24E04*(F0E)**2 

**  COMPUTE  MAX  DENSITY  OF  Fl-REGION  **«*•* 

_FNMAX*l .246C4* ( FOFl J *»2 

••  COMPUTE  MAX  DENSITY  OF  F2-REGI0N  •••*** 
tNK-A?=1.24E04*(  F0F2) **2 

**  COMPUTE  HEIGHT  OF  MAX  DENSITY  (SHIMAZ/.Kl  EQUATION) 

HMAX* 1490. 0/EM3000- 176.0 
**  COMPUTE  SCALE  HEIGHT  OF  F2-REGI0N  •****• 
TF*<AL0G(HMAX)/2.186E-02)-2C3.447 


3  •••• 


«*•* 
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c 

C******  HEIGHT  OF  F  1-REGION  SfcT  fQ  MIDPOINT  OF  E-  AND  F2-REGIONS 

C_  _  _  _  _ _ 

"  HMAXl=(HMAX*l 26.  t /2 . 

C 

C******  COMPUTE  SCALE  HEIGHT  OF  F1-R6GICN 
C 

TFl=(AL0G(HMAXl)/2.186t-02 1-203.447 
IF( IOPT.LT. 31  WRITE  (o,990)  I VB , I VE , l YR , IHR 
- 'Frit JPT.'E0.3i  liRlTE  1 6fT9T  YYRVmnThIYI  05",  1  hr  " 


C 

(>•••••  OUTPUT  SECTION  ••••» 

C 

C  " 

C*i«Mt  WRITE  HEAC^NG  FOR  SUMMARJT  PAGE _ •*•***  _ 

c  . . . .  ” 

990  FORMAT  I 1H1 , 26H IONOSPHERIC  PROFILE  VALID  , A6,2X,A6, 2X, I  2 , 1  5, 1HZ ) 
99  FORMAT  (IH1,24HI0NCSPHERIC  PROFTlF  FOR  t 3 1 3*1 5, IH2 1 


WRITE  (6,100)  CL  AT ,  NORS ,  CLONG ,  I  HEM ,  ( NAME  ( I  I  •  I  *1_»  4 ) 

100  FORMAT  flHO, 17HSTA1 ION  LOCATION'  ‘,'2TF6 •  1 ,  A  i  )  /  / 1 5X ,  4A6 ) 

I F < IOPT.EQ.l)  WRITE  (6,1000) 

"  ~  lRI0yT.~EQ.~2 )'  WRITE  (6,TgOn  _  '  .  . 

1001  FURMAT  ( lHOt^THTHIS  PROF ILE  BASED  UPON  LONG-TERM  DATA  TAPE) 

1000  FORM AT  (1H0,VX,  40HTHIS  PROF ILE  BASED  UPON  ITS  COEFFICIENTS) 

WRITE ~(6, 101 )  F0F2 , EM3000* FOE, F0F1 

101  FORMAT  ( 1  HO , 7HF0F2  =  , F S. 2 , 1 OX, 8HM3000  =  , F 5 . 2 , IPX , 6HF0E  *  ,F7.2, 
1 10X,  7HF0F  l  *  ,«=7.2» 

WRITE  (6,102)  GLAT ,  Z ANG,R 

102  FORMAT  ( 1H0, 23HGEOMAGNET IC  LATITUDE  *  ,F7.2,5X, 

121HS0LAP  ZENITH  ANGLE  =  ,F7.2//1X, 17H SUNS POT  NUMBER  *  ,F5.0) 
"WRITE  (6,103)  "TE, HE, ENE 

103  FORMAT  ( 1H0, 19HVALUES  FOR  6-REGION// 

1 5X, 1 5H SCAL  E  HEIGHT  *,F7.2,3H  KM/  f 
2 5X , 9HHE I GHT  =  ,F7.2,3H  KM/ 

35  X,  l'OHDEN  S  I  TY*  ,'FB  .0 , 13H  E  LECTh'okS/CC ) 

WRITE  (6,1040)  TF1 , HMAX 1 ,FNMAX 
1040  FURMAT  ( 1H0,20HVALUE S  FOR  F 1-REG  ION//  . 

1 3X , 15HSCALE  HEIGHT  =,F7.2,3HKM/ 

25X,9HHEIGHT  *  ,F7.2,3H  KM/ 

35X, 10HDENS ITY  =  ,F8.0,13H  FLECTRONS/CC ) 

- -p j n- T 6-.  1Q4,'  tf tflMA y , enma y - 

104  FURMAT  ( 1H0, 20HVALUE S  FOR  F2-KEGI0N// 
i$X,  15HSC4L1:  "HEIGHT  =,F7.2,3H  KM/ 

26X.9HHE IGHT  *  ,F7.2,3H  KM/ 

3  5  X ,  TO  HO  E  N  S  r  T  Y  *  "TTf.  O',  13h  EL  ECtROFfT/m 
£**««**  WRITE  HEAUING  FOP  PROFILE  PAGE  ****** 

- RirrmrsTTm - 

20 1  FORMAT  ( 1K1 , 1  OX , 24HELECTR0N  DENSITY  PROF  I LE//5X ,2HKM , 5X , 

nTHF-Trer,'  rGsrrsx  r«j  htt-  ptst  on  ,  5x  ,  shf  2-u  egi  w,TX7TH'Tcrr7cr,~7X', - 

210HCUMULATIVb,5X, I6HPLASMA  FREQUENCY, 5 X,5HSC ALE ) 

¥NP1=0,'  '  '  ~  -  - 

ENS AV=0* 


COMPUTE  AND  PRINT  VALUES  FOR  EACH  6  K.M  LEVEL 


C 
C 
C 

‘  Vi*95>.  ‘  ' 

LO  2C0  I-l«182 
J*l/2 

F2E=IH-HMAX1 1/TF1 
EZE*(H-HE  j/TE 
C 


£*■*•'* "'ELECTRON  0 E N is  I  T  I  ES  COMPUTED  "FOR  F2-RtGION  SISEO  "ON  A 
Cttuii  CONSTANT  SCALE  HEIGHT  IF  RELOW  THE  F2-PEAK  AND 

€••••••  UN  A  VARIABLE  SCALE' HEIGHT  IF  ABOVE  THE  F2-PEAK 

C 


’  IFlH.GT.HMAX)  Z EE  =  (H-HMAX ) /SCALE II) 

IF(H.LE.HMAX)  ZEE* ( H-HMA X ) / TF 

""  £E*ENE  *  EXP (  0 • 5* ( 1 VO-EZ  E-EXP  I -1 • 0*EZE I ) ) . 

EN*tNMAX*EXD  I  1 .O-ZE E-EXPI - 1 . 0*Z EE )) 

F  N = F  N  M  A  X  *  E  X  P  U .C-FZE-EXP (-1.0«FZt» ) 

ENP-EN+EE+FN 

IFlH.GT.HE.AND.H.LT.HMAXl.AND.ENP.LT.ENe)  CNP=EN£ 

IF (H.GT.HMAX 1. AND. H.LT. HMAX.AND.ENP.LT. FNMAX I  ENP*FNMAX 

I  F  (  ENP  7lT.  ENP  1.  AND .~H .  L  T  .'HMAX')  ENP  =  ENP1  . 

ENP1*ENP 

PLAS=8.97E-G3»SQRT( ENP ) 

STOKI I )*PLAS 

I F  ( MOD  (  1 , 2  )  .  EQ.T)  I  VAL  (  j )  ®PL  AS 
ENSAV*ENSAV*ENP»5.0EC9 

"  IF ( I . EQ. 1 )  GO  TO  200 . 

I  H*H 

WR  I T I  (6 , 202  )  IH  (  E  E",  FN,  EN ,  ENP ,  ENS  A  V ,  PL  AS ,  SC  AL  E  f  I ) 

202  FORMAT  ( 3X , l A, IX , 4< AX , F9.0 ) , 5X , iPE 16. 8 ,6X , 0PF7. 2, 7X ,F7.2  ) 
200  'H*H+8.0 

c 


C  *  * »  » * *  CALL  PLOt  “routine  If  REQUESTED  ~ 
c 

■" . lF(IPLOr.6T.O)  CALL  PLCTCvAU 

c 

. €***•*•  CALL  PUNCH  ROUTINE  IF  "RCQUbSTtf  •****• 

c 

. .  '  rrn pncr .xnvor  call  s  nrmrmvtNur,  rv'B ,  ivevtvrvtrri  ' 

c 

IF  IOPT  EQUALS '3,  ONLY  ONE  STATION  ANALYZED  PER  FUF2  AND  MI  3000 

.  I  F  I  I OPT. EC. 3  )  GO  TO  999 

C 

“  TTTWTT.-  ~i  N CRtN ENT  'RUTj'R — »ir.T. - - -  - 

C 

I HR* I HR* INC 

IF! IHR.GT.IEHRJ  GG  TO  999 
TC  8G0 
C 

-•  C»*  •  *v*  ~t X  I  T  "  AND  "TERN TNSTF  ION  . . — 


C 

777  CONTINUE 
ENCF1LE  9 
REWIND  9 
STCP 
end 
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tIBFTC  CURVE- 

SUBrGUT  INE  PLOT(VAL) 

INTEGER  GRIU(  46,91 ) ,  AHSC (46)  t  FREQ (46) 

"'INTEGER  EYE,BLANKT60I,UASH,fcV 
REAL  VAL ( 91) 

LATA  EYE, BLANK, DOT, LASH, EX/lHI,  li(  ,  1H.  ,  1H- ,  1HX/ 

L A  T  A ( ABSC ( I ), I  =  1 , 46) ✓ IH  , 1H  , 1H  , 1H  , 1H  ,1H  ,1HP,1H  ,1HL,1H  , 


1 ;HA , 1H  ,1HS,IH  , 1HM 

, IH  , 1HA, IH 

, IH  , IH 

, 1HF , IH  , IHR , IH 

,  1HE , 

2 IH  , 1HQ, 1H  , 1HU, IN 

, 1HE  j IH  , 1HN 

»1H  , 1HC 

,1H  , 1HY , IH  ,IH 

,1H  , IH 

)1H  , IH  , IH  ,1H  ,1H 

/ 

CATA(FKEQ( I ), I=l,46)/4H  0.0, 4H 

,4H 

,4H  1.0, 4H 

,  4H 

14H  2.C.4H  ,4H 

,4H  3.0, 4H 

,4H 

,4H  4.0, 4H 

,4H 

24H  5.0, 4H  ,4H 

, 4H  6.0, 4H 

,4H 

,4H  7.0,4H 

,4H 

34H  8.0,4H  ,4H 

,4H  9.0, 4H 

,4H 

,4H10.0,4H 

,4H 

44HI 1 .0 ,4H  ,4H 

, 4H12 ,0»4H 

,4H 

,4H13.0,4H 

,4H 

54H14.0,4H  ,4H 

, 4H 1 5 • 0/ 

CO  2GO  1*1,46 
CO  200  J=i,9i 
GRIL( I, J)*BLANK 
IFU.EC.l)  GKIO(T,J)*DOT 

IF1I.EQ.1. ANU. MOO (J,5).EQ»1)  GR ID ( 1 , J ) =£YE 
I  F  (  j" EC .1 6  f  GRI  6 1  I  *  J )  *DOT  ' 

I F  I  I  .  tQ .  16 . AND. MOU^J ,  5 ) . EQ. 1 )  GRID!  I  *  J ) =E YE 
"lFVi.EQ.3ij  GRIOII,  J)=OOT 

1 F ( I .EC.31.ANU.M0Dl J.5J.EQ.1)  GR  I D  l  l  ,  J  )  =E  YF 
I F ( I  *  EC .46 )  GRIOII,  J ) =00T 

I F ( I .EG. 46. AND. MOO ( J , 5 ) . EQ. 1 )  GR I D l I , J ) =E YE 

200  "CONTINUE  . . 

CO  300  1=1, 46 
'  GRIOII, 1)=EYE  ~ 

GRID! 1,91 )*EYE 
300  CONTINUE 

CO  1C  L=lj.91  _  _ _ _ 

" L 0 C  =  3 . C * V A Lf  j  ) ♦  1 . 6" ‘ 

10  GRIUILOC, IJ=EX  _ 

WRITE  (6,50 

5G  FORMAT  I 1H1.40X, lOHHt IGHT (KM) ) 

WRITE  (6,51) 

51  FORMAT ( 1HC, 11X, 3H10G, 2X , 3H 1 5C , 2 X , 3H200 , 2X , 3H2 50 ,2X, 3H300 , ' X, 

1  ^ 2 x~3H45 o ,  2 X ,  3H 4  5 0 ,’ 2  X  J  3 H 50 0  T 2 X ,  3 H 5 5 0,'2 X ,  3 H 6 0 Q ,  2  X ,  3H 6 5 0 , 2  X , 

2  3H700.2X,  3H750,2X,3HB00  2X,3H85C,2X,3H900,2X,3H950,2X,4H100C) 

CO  49  1=1,46 

WRITE  (6,30)  AHSC(I),FREQ(I),(GRID(I,J),J=1«91) 

SC  FORMAT  ( 5X , A 1 , 2X , A4 , 1 X , 91A1 )  . 

49  _  CONTINUE  _ _ 

RETURN  *  "  . 

END 
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ilBFTC  S1GNIF 

SUBXOUT I NE  SIGISTOR, IHGT. IVB, I VE. IYK, IHR) 
DIMENSION  S_TGR(182J,IHG_T(182),PLASQ(17) 
INTEGER  HE  IT ( 17) 

HE  I M 1 1  =  IHGT ( 2 ) 

HE  I T ( 2 )  =  IHGT ( 6  J 
HE  IT ( 15 )  =  IHGT ( 182) 

HE  II ( 16 )  =  IHGT ( 22  J 
P  L  A  SjQ_l  16)=STOK(22) 

--  Hfc  l  1 1  17  )  =  I HGT  C  32 )~ . * 

PL  ASQ( 1)  =  STOR ( 2 ) 

PLASQI 17 ) =STOR ( 32 ) 

PLASQ(2)=STGR(6) 

PLASQ(15l=STORri82r 
CO  I  I  * 1 f 181 

- mTrtmTrnEr.^TBRTi+i  nr  liirTo  v . . 

1  CONTINUE 

H  E 1 7  (  3 )  =  1 HGT ( 12) 

PL  ASQ  (  3  )_=  ST  OR  (  12) 

. k*i2"  . . . . 

GO  TO  5 

"4  HElTT3)*IHGTlT)  "  ' 

PLASQI 3)=STOR( I ) 

K*  I'  " 

5  UO _6  I=K,181 

TF ( ST ORff vn ; NE.S TOR (kn  GO  TO  7 

6  CONTINUE 

- - - 

PL  AS Q  (  4J  =  STOR  I  22) 

'  ~  K»22 

GO  TO  8 

7  he  1 1  ( 4  T  *  i  hot  rry . . 

PLASQI 4 )  =  STOR( 1 ) 

K=  I  . .  . . 

8  DO  9  I=K. 181 

'  . .  I F ( STOR! I ) .  LT. STOR ( I + 1)  )  GO  TO  9 

HE  I  T ( 7 ) * IHGT 1 1 ) 

. p La  s q (7 )= s  Tor (i )  '  . 

GO  10  10 

9  CONTINUE 

10  K=  I  -  5 

. . L*  I -10 

L  L  *  I  ♦  5 

”  HEIT(8)=THCT(LL)  '  ’  ’ 

HE  I  T  C  6 )  =  IHGT ( K ) 

- TTFITT5  7*1  HGTTU - - - 

PLASQI8)=ST0R(LL) 

PLASQI 6 )=S TOR l«J 
PLASQI 5 ) =STOR( L ) 

-  HElTl9)-lHGTT52r  - 
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PL  ASQ I  9 )  =  STOk ( 52 ) 

K  =  6  2 

IF(HEIT(9).GI .370)  K=72 

HE  I  T  ( 10  )  •'!  HGT  (  K )  "  . 

PLASU  10)  =  STOR(K) 

HE  I  T  ( 1 1 )  =  IHG T ( 62  ) 

PL ASQ ( II I -STOK ( 82  ) 

HEIT ( 12 )  =  IHGT ( 102) 

_  PLASQU2)  =ST0R(  10JJ _ _  _  _  __ 

HE  I r ( 1 3~f=  I HG f ( 132 ) 

PL  A  SO  (  1 3  )  *STOP.  C 132) 

HtIT(14)=IHGT<162) 

PLASQ(14)=ST0R(162) 

15  GO  100  1=2,17 

I F  ( Hb  1  T I  I  l.GE.HEIK  1-1) )  GO  TO  100 
‘101  I T  EMP  =  Htl T (  I  ) 

HEII  m-hE!Tt  I- 1 1 
HE  I T ( I - 1 )  =  1  TEMP 
TEMP=PLA$Q( I  ) 

PLASQ ( I ) -PLASQ ( l-l) 

PLASQ  (  1^1 )  *TfcMP _ _ _  _  _ 

. CO  TO  15 . 

100  CONTINUE 

IHR1=  IFR-70 

IFUHK1.LT.0)  IHR1=2330  _ 

IHR2MHROO 

_IHR_l=JHKl  +  100gC _ _ __  _ 

WRITE  "(9,991  IVB,  IYR,  l  VE*  IYR*  IHR1 »  I  HR  2 
WRITE  (6,98)  IVB,  IYR,  IVE_*  £YR,  IHRI, IHR2 
99  FORMAT  (6HVAL11)  , A6*13,3H  -  ", A6» 1 3, IX, l4, 5HZ  TO  ,14, 

1 32hZ  REMOVE  THIS  CARD  BEFORE  USING) 

98  FORMAT'  ( 1H1 , 26H IONQSPHER I C"PROF ILE” VAL I D  ,A6,I3,3H  -  ,A6, 13, IX, 14, 
15HZ  TO  »J4_,1HZ)  _  __ 

""GO  20 0  1  =  1,17" 

JWR  I  T E  (  6j  9 7  )  HE  I  T  (  IJ  ,  PLASQt  I  ) 

WRITE  (9,9*6)  HE  I T  (  I )  *  PLASQ  ( if 
200  _  CONTINUE 

97  FORMAT  *  (  lX,4HICNCVi2X,  1*4,  F12. 2) 

96  FORMAT  ( 4H 10N0 , 12 X , 1 4 , F 1 2 . 2 ) 

"'"RETURN*”"  "  . 

ENG 


Appendix  B 


Sample  Computer  Output 
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UNCLASSIFIED 


IONOSPHERIC  PKPEILL  VALID  23  HAY  07  JUN  70  19007 _ 

STATION  LUC. AT  I  ON  2  4.  ON  RS.Ow _ 

_  EGL 1  N  FANGE _ 

THIS  PROF  11  E  BASED  UPON  ITS  COEFFICIENTS _ 

F0F2  °  3,25 _ M3000  *  2.76 _ FQF  =  4.04 _ F0F1  _= _ 5.59 

GEO.UO.NfcTlC  EAT1TU0E  =  35.48  SOLAR  ZENITH  ANGLE  *  18.11 _ _ 

SUNSPOT  NUMBER  =  SO. _ _ 

VALUES  FOR  F  -REG  ION _ 

SCALE  HEIGHT  »  15. 56  KM 

HFTiHT  »  120."CR)  K*  '  '  "  . *'  . 

DENSITY  *  202076.  ELECTRUNS/CC _  _  _ 

VALUES  FOk  Fl-MEGION _ 

SCAl.  £  HE  IGHT  =  47.  56  KM 

HEIGHT  *  '24U34'kT7  '  *  "  . . 

OENSITY  « _ 3B6SBS.  ELECTRONS/CC _ _ _ 

VALUES  FUR  F  2-RF.G  1  OlM _ 

SCALE  HEIGHT  *  66.21  KM 

- HTTCTHr. — 36XFrK*r - — - — -  •— 

OENSITY  =  106182“ ,  ELECTRONS/CC 
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